function [v,d] = getranddraw(k2,income, ns)
% prepare antithettic draws
rng(175);
ns2=ns/2;
v1=randn(ns2,k2);   
vminus=-v1;
v=[v1' vminus'];
d=randsample(income,ns)';
end